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ABSTRACT 

We study the mass distribution in the late-type dwarf galaxy NGC 2976 through stellar kinematics obtained 
with the VIRUS-P integral-field spectrograph and anisotropic Jeans models as a test of cosmological simu- 
lations and baryonic processes that putatively alter small-scale structure. Previous measurements of the Ha 
emission-line kinematics have determined that the dark matter halo of NGC 2976 is most consistent with a 
cored density profile. We find that the stellar kinematics are best fit with a cuspy halo. Cored dark matter halo 
fits are only consistent with the stellar kinematics if the stellar mass-to-light ratio is significantly larger than 
that derived from stellar population synthesis, while the best-fitting cuspy model has no such conflict. The 
inferred mass distribution from a harmonic decomposition of the gaseous kinematics is inconsistent with that 
of the stellar kinematics. This difference is likely due to the gas disk not meeting the assumptions that underlie 
the analysis such as no pressure support, a constant kinematic axis, and planar orbits. By relaxing some of 
these assumptions, in particular the form of the kinematic axis with radius, the gas-derived solution can be 
made consistent with the stellar kinematic models. A strong kinematic twist in the gas of NGC 2976's center 
suggests caution, and we advance the mass model based on the stellar kinematics as more reliable. The analysis 
of this first galaxy shows promising evidence that dark matter halos in late-type dwarfs may in fact be more 
consistent with cuspy dark matter distributions than earlier work has claimed. 

Subject headings: galaxies: individual (NGC 2976) — galaxies: dwarf — galaxies: kinematics and dynamics 
— galaxies: spiral — dark matter 



1. INTRODUCTION 

The observations of kinematics in low surface bright- 
ness (LSB) and dwarf-late type galaxies have stubbornly re- 
sisted giving clear evidence for the cuspy Navarro-Frenk- 
White (NFW) dark matter (DM) halo pr ofiles that N- 
body s imulations with ACDM inputs predict (iNavarro et alJ 
1996b). Instead, most LSBs and late type dwarfs suggest 
cored DM halos (e.g. lOh et al.1 120081: Ide Blok et alj 120081 : 
iKuzio de Narav et alj2008tlSpano et al.l2008tlOh et alj201oi) 
or the observations are not constraining eno ugh to rule out 
cusps (ISwaters et alJl2003l ISimon et alJl2005l) . Some simula- 
tions have produced cored DM halos by rapidly removing the 
baryonic dis k which causes the D M halo to expand to a cored 
equilibrium (Navarr o et al.ll996al) . initializing numerical sim- 
ulations with a pri mordial bar that forms a resonance with and 
disrupts the cusp (Weinberg & Katz 2002[), or by implement- 
ing a supernova feed back recipe in high-re solution hydrody- 
namical simulations (Governato et al. 2010). This puzzle has 
also motivated proposals for additional particle properties of 
dark matter beyond the weakly intera cting, cold paradigm 
such as collisional dark matter (e.g. Spergel & Steinhardt 
2000), warm dark matter (e.g. Hog an & D alcanton 2QQ0|) , and 
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ultra- l ight boson fluid models ("e.g. [P eebles! [20001 : iGoodmaiil 
I2000t iRindler-Daller & Shapirdl201 ll) . Other questions crit- 
ically rely on knowing the DM density structure in galax- 
ies, such as the prospects of DM annihilation searches that 
depe nd sensitively on the true density profile in galaxies 
(e.g. iDiemand et aT] 120081) . Most of the extant attempts to 
determine DM radial profiles rely on gas as the dynami- 
cal tracer. A number of works have studied nearby disk 
galaxy kinematics with longslit stellar kinematics to infer 
DM h alo profiles jCorsini et alj|1999t ICorbelli & Walterbosl 
120071) or stellar mass-to-light ratios for isothermal sheet 
mode l s (Ivan der Kruit & Freemanll 984: Bah call & Casertanol 
| 1984t Ivan der Kruit & Freemanl ll986riBottema et all 1 1987b 
Bottemall 1 989bllaL 1 1 9901: iBottema et alJl!991b lBottemalll99a 
ISwaters! 119991) . but bett er structure constr a ints come from 
2D s p ectroscopy (e.g. ICopin etall l2004t iKrainovic et alJ 
1 20051: ICappellariet all 120061: Ivan den Bosch et all 120081: 
Weijmans et alJl2009tlMurphv et alj|201 ll) . 

NGC 2976 has made one of the c leanest cases for a cored 
DM halo via its gaseous kinematics (ISimon et al.ll2003l here- 
after SBLB03). In our first attempt to derive DM mass pro- 
files from stellar kinematics we chose NGC 2976 due to 
several of its properties. (1) NGC 2976 is an SAc dwarf 
galaxy in the M81 group. There are some dynamical in- 
dications ( Spek kens & Sellwo od 2007) and photometric in- 
dications (Menendez-Delmestre et al . 2007) that a weak bar 
may be pres ent, but NGC 2976 is usual ly given an unbarred 
designation dde Vaucouleurs et al][l991l) . (2) NGC 2976 has 
some dark patches that are likely due to dust, but its dust 
content is modest for its Hubble class and distributed rather 
regularly. A full treatment regarding the potential impact of 
dust on the measured kinematics is beyond the scope of this 
work, but several liter ature estimates of the dust content exist. 
IWilhams et al.1 (120101) fit star-formation history models from 
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the broadband colors of resolved stars in the Advanced Cam- 
era for Surveys Near by Galaxy Survey Treasury (ANGST) 
(Dalcantonet al. 2009]) by modelling 0.8 magnitudes of dif- 
ferential extinction in the V-band for ages above 100 Myr, 0.5 
magnitudes for younger ages, and a foreground screen of 0.46 
magni tudes. Our data lie w ithin the "INNER- 1" region of that 
work. iPrescott et al.l (120071) use Spitzer 24-fim data to estimate 
As ~1.5. The lowest estimates comes from S BLB03 with 
Ag~0.23 based on an inclination proscription (Sak ai et ail 
2000). Although there is a large range in the estimated ex- 
tinctions, the values from each method are on the low end for 
the late-type dwarf population. (3) NGC 2976 appears to be 
dark m atter dominated at r >500 pc according to SBLB03 (al- 
though lde Blok et ail (120081) disagree), so the impact of stellar 
population synthesis (SPS) mass-to-light (T„) uncertainties is 
minimized in this target. (4) The stellar surface brightness 
profiles in many bands are very smooth and indicate a bulge- 
less disk with a small nuclear star cluster and a break to an 
outer disk at -1.2 kpc (SBLB03). 

The distance to NGC 2976 fro m the tip of the red gian t 
branch method is 3.56±0.38 Mpc (iKarachentsev et al.l 2002); 
we adopt a distance of 3.45 Mpc and a scale conversion of 
16.7 pc arcsec -1 in this work for consistency with SBLB03. 
The total mass of NGC 2976 is estimated to be 3.5 xlO 9 
Mq based on the inclination-corrected line width of 165 km 
s _1 ( SBLB03). The inclination is variously estimated as 
61?5 (Ide Vaucouleurs et al.| [l99l. 61?4 (SBLB03), and 64?5 
(Ide Blok et al.ll2008h. TheHI heliocentric velocity is 4.0±2.0 
km s _1 dltil & Israel 2002b). The stellar velocity dispersion 
in NGC 2976 has not previo usly been repo rted. Measure- 
ment attempts were made in iHo et ail (|2009), but the results 
were unresolved in the presence of the best <Ji ns t — 42 km 
s _1 instrumental resolution and unreported. They estimate 
a =36.0±16.8 km s _1 from a correlation for their sample 
with the measured [NII]6583A line width. 

We present observations of NGC 2976 with the large field- 
of-view fiber fed Visible Integra l field Replicable Unit Spec- 
trograph Prototype (VIRUS-P) (iHill et all 12008) to concur- 
rently measure the gaseous and stellar kinematics, fit mass 
models, and study the dark matter halo profile shape in 
the context of the "core-cusp" controversy with a collision- 
less tracer. Our data reduction and kinematic measurements 
are described in $2\ We fit the stellar kinematic data with 
anisotropic Jeans models in $4] We perform fits to our [Oil] 
rotation curve and the SBLB03 Ha rotation curve in |5] In Sj3] 
we investigate the constraints on T* through matching SPS 
models to our spectral data and its effect on the mass models. 
Finally, we review our conclusions in $6] As is customary, all 
values of T* are given in solar units for the indicated band 
with the "*" indicating a mass exclusively for the stellar com- 
ponent. We adopt an absolute solar mag nitude in the R-band 
ofM_0 u = 4.42 throughout this work (Bin nev & Me rrifield 
[19981) .' 

2. OBSERVATIONS AND DATA REDUCTION 

Over April 27 through May 1 of 2009, we took 18 hours 
of science pointings on NGC 2976 with the VIRUS-P in- 
strument and the 2400 lines mm -1 VP2 grating on the Mc- 
Donald Observatory's 2.7m Harlan J. Smith telescope under 
non-photometric conditions. Transparency was continuously 
monitored by stars in the guider camera's data, which was 
read out and saved every five seconds and used to make rel- 
ative flux calibrations. The seeing ranged from l'.'4-3'/0 full- 
width-half-maximum (FWHM). The VIRUS-P field-of-view 



covers 1.'61 x 1 f 65 with fibers projected to on-sky diameters 
of 4'.'235. VIRUS-P has a one-third filling factor, so we spread 
three dithers across the galaxy to maintain continuous spatial 
coverage in the middle and to maximize the spatial extent. 
We set the instrument to measure 3680-4400A with R=2400, 
nominally. The spectral data were taken under lxl binning, 
yielding four pixels sampling the instrumental FWHM. In- 
dividual exposures were of 30 minute duration. We took 
20 minute sky nods offset from the galaxy by 10' between 
the science data frames. Our best stellar kinematics come 
from around the G-band at 4300A and its many surround- 
ing, mostly Fe, features. We also measured the [OII]AA3726, 
3729 lines and all Balmer lines higher than and including H7. 
T he data were reduce d with the vaccine pipeline described 
in lAdams et al.l (1201 ll) . We have corrected all observations 
to the heliocentric frame. The instrumental wavelength zero- 
point is observed to drift by —10 km s _1 over normal swings 
in nightly operating conditions. We track this and correct the 
zeropoint by fitting the 4358. 3A Hg I skyline in every nod- 
ded sky exposure. Quoted wavelengths are not corrected to 
vacuum conditions. 

2.1. Binning 

We have, in total, 738 spectra with signal-to-noise (S/N) 
ranging up to 60 per pixel. In order to extract reliable stellar 
velocity dispersions under high S/N conditions, we bin our 
data as shown in Figure Q] Constant sized bins were chosen 
for simplicity. The stacks are made with one interpolation 
to a common, linear wavelength scale considering the differ- 
ent wavelength solutions in each fiber and the individual he- 
liocentric corrections. Membership in a bin was determined 
solely from each fibers' central position; fibers extended over 
multiple bins were not given partial weights. We tested ex- 
tractions using partial weights, and the results agreed to those 
presented to within the estimated errors. This is expected 
since the square bins are considerably larger than the fibers, 
and the usual fractional area that seeps in and out of the square 
bin for a collection of fibers is only 10%. The spectra were 
not degraded to have a common instrumental resolution prior 
to binning. We made this choice to minimize the covariance 
between spectral channels. Instead, the average instrumental 
resolution was associated to each bin. We compared the kine- 
matics derived from making resolution-matched stacks and 
from simply degrading the kinematic template spectra to the 
average instrumental resolution in each bin. The two meth- 
ods agreed, on average, and only displayed differences within 
the estimated errors. Adaptive binning strategies were not 
adopted, so there is a significant range of S/N in the binned 
spectra. The mean S/N per pixel is 29 with a dispersion of 23. 
We test wh ethe r this range in S/N influences our extracted ve- 
locities in £|2.4l and find no impact. In the case of NGC 2976, 
adaptive spatial binning is not ideal and would lead to the loss 
of spatial resolution in the outer parts where the contribution 
of Vios to V rms is significant. The analysis of the stellar kine- 
matics is made with these 16" x 16"bins, while the analysis of 
the gas kinematics is made under the native fiber diameters of 
4'.'235. 

2.2. Spectral resolution 

Our observations were taken with an instrumental disper- 
sion of between 40-60 km/s across different fibers and wave- 
lengths. We measure this resolution to 0.5 km s _1 uncertainty 
in every fiber and wavelength with our twilight flats and high 
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resolution solar spectrum data dKurucz et alJI 1984b . Finally, 
the binning process combines fibers with different instrumen- 
tal resolutions, per our dithering scheme. We average the in- 
strumental resolution for each bin. A fit to one fiber is shown 
in Figure [2] The fit to all bins is summarized in Figure [3] 
The solutions agree with arc lamp data, although the sparse 
number of available arc lamp lines leads to a less constrained 
solution. We also took spec tra of a number of template stars 
(Prugniel & Soubiran 2001) to test our instrumental resolu- 
tion and find agreement. We experiment with degrading all 
the data in a bin to the maximum instrumental resolution prior 
to stacking, but we find differences only at levels far smaller 
than the formal errors. 

2.3. Gaseous kinematics 

We measure the gaseous kinematics through the 
[OII]AA3726.032, 3728.815 doublet. Without binning 
and in each fiber, we simultaneously fit two Gaussian func- 
tions over an 18A window. Five parameters are fit through 
a least-squares minization: the intensities of each emission 
line, a constant continuum, the radial velocity, and the line 
width. The best fit models are perturbed with the estimated 
flux uncertainties in Monte Carlo realizations to generate 
velocity uncertainties. The median measured intrinsic line 
width is a = 20.4 km s _1 with no strong spatial gradients. 
We make models of the circular velocity profile from the 
line-of-sight velocity measurements in |5] 

2.4. Stellar kinematics extraction 

We fit the stellar kinematics in each bin with a maximum 
penalized likelihood estimate of the Gaussian line-of-sight ve- 
lo city distribution (LOSVD) in pixel space via code described 
in lGebhardt etafl (12000h. We use the empiric al, R = 10k stel- 
lar r Template7oflPn^ (ELODIEv3.1) 
convolved to the instrumental resolution of each bin. All data 
and templates are normalized prior to convolution by running 
two boxcars over each spectra. The first boxcar has a 40A 
width. All pixels deviating from the smoothed continuum es- 
timate by > 1.5cr are masked in the second boxcar, which has 
a 14A width. The specifics of this normalization have little 
impact, as judged by varying the normalization parameters, 
on the convolution since both the flux-calibrated VIRUS-P 
and ELODIE spectra hava very shallow slopes over the ex- 
tracted wavelengths. 

The choice of templates is made from amongst the 1959 
available, although only a small number are required to de- 
scribe the data. The membership is determined by manually 
iterating the list to find a local miminum in root-mean-square 
(rms). Several such minima can be found, but the exact choice 
is unimportant so long as similar stellar types are included. 
We list our chosen, final set in Table Q] A Wolf-Rayet star 
takes on a significant weight. However, its only function over 
this wavelength range is to lower the combined equivalent 
width of absorption features and is degenerate with, say, the 
A7m star HD003883. The metallicity has been modeled by 
ANGST at [M/H]=-0.12 for this region and by imposing a 
constraint that metallicity grows with time. Most of the F and 
G giants in our template list have higher metallicity. How- 
ever, experiments including the template HD 148856, a G8III 
at [Fe/H]=-0.26, do not improve the fit. 

Although we track formal uncertainties, we also tally an 
empirically determined, systematic uncertainty as the rms 
from the best-fit model. The two generally agree, expect at 



high S/N where template mismatch becomes evident. All fur- 
ther uncertainties are based on the empirical uncertainty, al- 
though these may overestimate the actual uncertainty. Un- 
certainties in velocity are determined by making Monte Carlo 
realizations of the best-fitting model with simulated noise de- 
termined from the residuals. 

Tests are run to determine the limits of reliable kinematic 
extraction under a range of S/N and intrinsic velocity disper- 
sions. First, a subset of the template spectra are combined, 
then convolved by a simulation dispersion and the instrumen- 
tal resolution, and noise is added. The extraction of the veloc- 
ity dispersion is shown in Figure|4] The errors are accurately 
estimated, and no systematic effects are seen down to S/N>5 
and a > 10 km s _1 . Next, a similar test is made to capture the 
possibility of template mismatch. We combine the ELODIE 
templates for HD000432, HD068380, and HD081809 in a 
21:53:27% ratio. The stars are not part of our fitting tem- 
plate, but they have similar spectral types. The extractions of 
velocity dispersion are shown again in Figure [4] The errors 
are marginally larger, but again they are accurately estimated 
and without systematic trends. 

Representative spectra and their best-fit models are given 
in Figure [5] Several additional corrections are made. The 
instrumental resolution uncertainty, estimated as 0.5 km s~\ 
is propogated as a random error to V rms along with the sta- 
tistical errors. The ELODIE headers quote broadening of the 
stellar features for some stars beyond their R = 10k resolution; 
we are not certain whether this broadening is physical, such 
as by binaries, or a spurious artifact. The average value for 
our templates is a — 5.2 km s~ x . We treat the broadening as 
physical and subtract off the ELODIE instrumental resolution 
and the average broadening in quadrature from our VIRUS-P 
instrumental resolution prior to template convolution. How- 
ever, the effect is small compared to the final uncertainties. 
Additionally, the systemic velocity is estimated directly from 
our data. The inverse-variance weighted average of our stellar 
bins is 4.60 km s _1 , con sistent with an earlie r optic al deter- 
mination of 3 ±5 km s -1 Ide Vaucouleurs et al.l d 19911) a nd the 
HI measurement of 4±2 km s _1 IStil & Israeli d2002bt) . This 
is subtracted from the line-of-sight velocities before forming 
V rms . The line-of-sight velocity and velocity dispersion for 
each bin, with uncertainties, are given in Table [2] Figure [6] 
shows the maps of line-of-sight velocity, velocity dispersion, 
and uncertainties that we measure from the data. The final 
velocity dispersions are not allowed to fall below 10 km s _1 
in order for their error estimates to impact the error on V rms . 
This is consistent with the lowest velocity dispersions that we 
can reliably extract in simulated data. 

Jeans models make predictions for the projected, second 
velocity moment. For a Gaussian kernel, the second-moment 
is simply V rms = \/V l I os + a 2 . The chosen spectral win- 
dow fit is 4110A< A <4340A which includes the G-band , 
a strong Ca absorption line at 4227A dWorthev et alJI 19941) . 
and a large number of weak Fe features. Fits to the spectral 
window 3980A< A <4100A give similar measurements, al- 
though noisier due to the smaller bandpass and less prominent 
features. We show the consistency between spectral regions 
in Figure [7] 

A careful viewing of the velocity dispersion map in Fig- 
ure [6] shows a profile that is nearly flat, but perhaps contains 
higher velocity dispersions in the northerly direction. This 
pattern remains whether we remeasure the instrumental reso- 
lution within each night's data or adopt a resolution fixed with 
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time. However, the evidence for a velocity dispersion gradient 
is within the noise. 

3. STELLAR POPULATION SYNTHESIS CONSTRAINTS 

A loose but independent estimate on T* can be made 
by comparing stellar population models to Spectral En- 
ergy Distribution (S EP) data, either photometrically (e.g. 
IB ell & de Jongll200ll) or spectrophotometrically. Substantial 
systematic uncertainties in, for instance, the initial mass func- 
tion (IMF) and the properties of therm ally pulsing asym p- 
totic giant branch (TP-AGB) stars (e.g. iMarastonl 120051) at 
NIR wavelengths limit the precision of M/L constraints from 
SED fits. Experiments designed to provide optimal obser- 
vational constraints on stellar M/L's by applying isothermal 
sheet relations are un derway dHerrmann & Ciardullol 120091: 
iBershady et alJ2010allbh . but the goal of this work is instead to 
model and fit all mass components simultaneously. Specific to 
NGC 2967, SBLB03 found tension in the maximal disk value 
of T*,k < 0-09io'og compared to the higher values implied 
by some SED fits (one of their models has T*. —2). 

We analyze the stellar population by fitting spec- 
tra with stellar p o pulatio n synthesis (SPS) models from 
iBruzual & Chariot! d2003l) and the preliminary release of 
their 2007 ver sion that incorporates ne w TP-AGB values 
dBruzua]||2007l). We u se the Padova 1 994 dAlongi etalJI 19931: 
Bress an etaTTl 1993b iFagotto et ail 1 1 994allbh iGirardi et all 
19961) and Marigo 2007 (Marigo & Girardill2007l) evolution^ 
ary tracks, respectively. Both Chabrier and Salpeter IMF's 
are tried, which represent reasonable lower and upper bounds 
on T,. To achieve robus t star formation histor ies, we use 
the same 39 templates as iTremonti et al.l ((2004) which en- 
tail combinations of three metallicities (Z — 0.2, 1,2.5Z Q ) 
and a variety of star formation histories (instantaneous bursts 
of age 0.005, 0.025, 0.10, 0.29, 0.64, 0.90, 1.4, 2.5, 5, and 
1 1 Gyr; a 6 Gyr old population under a constant star forma- 
tion rate, and two tau models sampled at a 12 Gyr age with 
tsfr = 5,9 G yr). We further add a grid of dust extinction 
with the form of lCalzetti etail (120001) over 21 values of E S (B- 
V) uniformly spaced from to 1. We resample our spectra 
to the rest frame, convolve the templates to match the 8A 
FWHM SINGS instrumental resolution, and mask out win- 
dows 4A wide around each of the Balmer lines, [Oil], and 
[NeIII]3869. Repeated observations of spectrophotometric 
standards stars with VIRUS -P have shown the relative flux 
calibr ation to be accurate to better than 10% (Ada ms et al.l 
l20Tlh . The match to the SINGS 20"x20"drift-scan spec- 
trum is excellent as shown in Figure [S] Due to the larger 
bandpass, we quote values by fitting to the SINGS spectrum. 
Similar, but less constrained, values result from the VIRUS- 
P data. The templates' normalizations are fit through least- 
square minimization, first individually and then in all 334,971 
two-component combinations from the metallicity, star for- 
mation history, and dust grids. The relative probab ility of each 
mode l is calculated as exp(— x 2 /2) (e.g. [kauffman n et al.1 
120031) . All quoted T* values include the effects of dust. The 
la confidence intervals centered on the highest probability 
T*. R are 0.63±0.39, 1.23±0.52, 1.04±0.23, and 1.42±0.42 
for the BC03/Chabrier, BC03/Salpeter, CB07/Chabrier, and 
CB07/Salpeter models respectively. Given the tight ranges 
from the statistical errors alone, the M/L uncertainty 
is dominated by systematic uncertainties. We consider 
^*.R =1-1 ±0.8, the mean and the symmetric uncertainty 
that encompasses the la confidence intervals for all four es- 
timates, as the best spectrophotometric limit; therefore we 



present mass models with T*,ij both freel y fit and penalized 
by this conservative confidence interval in i j4.3l andl5l 

Finally, we investigate the Y„ limits enabled by broad-band 
photometry. We use the optical-through-2MASS datapoints 
of SBLB03 by assum ing 10% errors and two IRAC datapoints 
of iDale et all (12007b for NGC 2976 without any aperture cor- 
rections as shown in Figure [9] The same st ellar population 
mode ls are fit through the EAZ Y package ( Brammer et al. 
2008). The best-fit value of Y„ is consistent with our spec- 
tral fits. Ho wever, it is starkly inconsistent with the SPS 
value u sed by lde Blok et al.l d2008) and calibrated in lOh et all 
(2008) that renders NGC 2976 to be dominated by the bary- 
onic mass. The relation derived there is based on SPS models 
that have aged 12 Gyr and does not match the observed colors 
of NGC 2976. Our fit to the S EP proscribes TV B = 0.64 
and T* i3 . 6iUm = 0.18 while the Ide Blok et all (12008) model 
proscribes T„^, 6fim = 0.66. 

4. JEANS MODELS 

The minimally necessary components to the NGC 2976 
mass model are a stellar disk with a spatially uniform mass- 
to-light ratio, atomic hydrogen, and a dark matter halo under 
a power-law parameterization. We use extant photometry to 
infer the distribution of the first two, and we use kinematic 
measurements and Jeans model solutions to infer the latter. 

4. 1 . HI deprojection 

We use the robust weighting (1Z = +0.5), zeroth moment 
map of 21cm atomic hydrogen i n NGC 2976 from The HI 
Nearby Galaxy Survey (THINGS. IWalter et alj|2008h to char- 
acterize the HI mass model. The HI distribution is highly 
clumped around the two off-center star forming complexes 
and unlikely to be axisymmetric. Nevertheless, w e perform 
Multi-Gaussian Expansion (MGE, Cappellari 2002) fits to the 
data, subject to the axisymmetric limitations of this work's 
modelling an d calculate a resu ltant circular velocity profile 
(Appendix A. ICappellarill2002l) . SBLB03 have shown that HI 
is a dynamically somewhat important component at r > 80", 
but that H2 (from CO measurements) is a minor contribu- 
tor to the gravitational potential at all radii. Therefore, we 
neglect the mole cular component. SBLB03 fit the HI from 
an older dataset dStil & Isr ael 2002a) and by assuming an in- 
finitely thin disk; they present a circular velocity profile that 
is in general a factor of two times larger than our deriv ed val- 
ues. We apply the MGE mass model to the fits of i!4.3l and Sj5j 
although the presence of the HI component does not strongly 
influence our final results. The contribution from HI to the 
circular velocity is given in Table [3] 

4.2. Stellar deprojection 

We use the R-band image taken at the Kitt Peak Na- 
tional Observatory's 2.1m telescop e from the Spitzer Infra red 
Nearby Galaxies Survey (SINGS. iKennicutt et "all 120031) to 
model the stellar mass distribution. The MGE fit is shown 
in Figure[l0]with the terms listed in Table|4] The MGE model 
fits both the nuclear star cluster and the inner and outer disks 
well. The inferred vertical-to-radial scale length is 1:8 at a 
nominal i = 63° over most of the radial range, although it 
is near-circular for the nuclear star cluster. Strong color gra- 
dients are not seen in NGC 2976 (SBLB03), so we limit our 
analysis to a single value of T^.^ for all components. There 
is some debate as to the best filter to use for accurate recovery 
of stellar mass. Thermally pulsing asymptotic giant branch 
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stars and polycyclic aromatic hydrocarbons disfavor the NIR. 
Dust extinction, varieties of star formation history, and nebu- 
lar emission dis favor the opti c al. Th e studies on NGC 2976 
of SBLB03 and lde Blok et al.1 (120081) usedK s and Spitzer 3.6 
/im data, respect ively. iPortinari et al.l (12004) advocates i-band 
photometry and Zibetti et al. (2009) calibrates for combina- 
tions of i and H band photometry. We chose to use the avail- 
able R-band image primarily for its depth and resolution in 
this work. However, that our fits to the spectral energy distri- 
bution are consistent across multiple bandpasses ([J3]l means 
this choice is unimportant for the analysis of NGC 2976. 

4.3. Best-fitting dark matter halo 

We use the Jeans A nisotropic MGE modelling package 
(JAM Cappellarl 120081) to fit the binned stellar kinematic 
field in NGC 2976 in lieu of more computationa lly inten- 
sive Schwarzschild modelling ( Schwarzschil dl 1 9791) . The fits 
are made to the projected, second-moment velocity (V rms = 
\/V 2 os + a -2 )- We make models assuming a single anisotropy 
parameter ((3 Z = 1 — {<j z /<Jr) 2 ) and a spatially constant Y*. 
The spherical DM halo's radial profile is approximated as a 
power law 

p(r) = po x (r/1 pc)- Q (1) 
while the NFW function has the form of 



p(r) 



(r/r a )(l + r/r s ) 2 



(2) 



with r s being a scale radius and pi being a density related to 
the critical density and the halo overdensity. 

The power-law approximiation is justified since the core ra- 
dius of NGC 2976 is likely to he at r s ~ 2.5 kpc (SBLB03, 
Appendix B), and our stellar data do not extend into the 
asymptotic region of the rotation curve. This choice also 
aids comparison to previous work whereas SBLB03 also used 
this parameterizatio n and most models for NGC 2976 in 
Ide Blok et al.l ([2008 ) minimized to solutions with lower lim- 
its on radial scale parameters that then reduce to the power- 
law form. We place the JAM code within a non-linear least- 
squares minimzation package (MPFITv . 1 . 643) to reach our 
optimal solutions. 

The best-fit, five-parameter model has a shallow minimum 
at T, >fl =3.49, i =63.3°, & =0.432, p (r = lpc) =0.260 
M Q pc" 3 , and a =0.235 at xl=87 =77.0. The second- 
moment model and residuals are shown in Figure QT| The 
enclosed mass distribution for this model is shown in Fig- 
ure [12] The anisotropy we measure is larger than the (3 Z ~ 
0.3 commonly fit by the same method i n E's and SO's 
(ICappellari et al.ll2007t IThomas et all 120091) and slightly be- 
low the J3 Z ~ 0.5 — 0.8 range commonly found in large 
spirals (Ivan der Kruit & de Griisl [19991; IShapiro et all 120031; 
iBershadv et al.ll201 ll) . We test an isotropic solution as well, 
and find a similar solution in the remaining parameters. The 
kinema tically determined inc l inatio n is marginally above the 
value in de Vau couleurs et al.l {l 99 1 ) based on photometric el- 
lipticity (i = 61.5°). However, we find i =64.6±1.5° from 
our [Oil] tilted ring [Oil] fit |5] so the kinematically deter- 
mined value is reasonable . Inclinations are ofte n poorly con- 
strained by kinematic fits (Kra inovic et al.l2005l) and can have 
a strong degeneracy with f3 z . There is tension in the T* value 
compared to that which we determine from the stellar popula- 
tion synthesis (SPS) fits presented in S|3] Using the constraint 

1 http: //cow .physics .wise . edu/~ craigm/idl/idl .html 



we derive from SPS, we add an additional x 2 term containing 
AT* fl=0.8 with a central value of T* v r=1.1 to the \ 2 values 
from the kinematics per 



x 2 = 



AT* 



E 



^i,rms ^i, model 



(3) 



The parameters then reach a minimum at Y* #=1.158, 
i =65.0°, B z =0.450, p (r = lpc) =45.7 M pc~ 3 , and 
a =0.90 at xt=ss =77.1. By marginalizing over the other 
four parameters, the data disfavor a=l at 0.8er significance, 
a=0.6 at lcr significance, and a=Q at 2.2er significance. The 
Vrms map and residuals are shown in Figure QT| The en- 
closed mass model is shown for the joint constraint in Figure 
[T2l This solution is within the la confidence interval from 
the kinematic data only. Finally, we fix T*./?=l and a =0.1 
as an illustration of a cored DM model. The V rms map and 
residuals are shown in Figure [TT] After optimizing the re- 
maining three variables, this model yields xt=88 =93.9 and 
is excluded with high confidence. The crucial difference be- 
tween this final model and the data is in the generally more 
circular V rms contour in the former. 

4.4. Parameter degeneracies 

We here investigate the degeneracy between a and T*. By 
our parameterization, there is an obvious degeneracy between 
po and a with a weaker degeneracy on T* (Figure [T3l). The 
most important degeneracy for the purpose of constraining the 
mass budget is that betweeen T* and a. In NGC 2976, they 
anti-correlate. A similar exercise is done with the joint kine- 
matic and SPS likelihoods by which the cored DM fit is ex- 
cluded at 2a significance but a pure cusp is permitted. We 
conclude that honoring even a loose T* limit makes a DM 
halo measurement entirely consistent with the NFW form. 
The cored model is disfavored at modest significance but 
amenable to stricter limits through more extended instrument 
pointings and higher S/N data. 

4.5. Other DM distributions 

We apply the power-law function as the primary DM den- 
sity distribution, but we briefly test and discuss alternatives. 
The pseudo-isothermal function is commonly applied to DM 
halo data as a well-fitting cored model, although it lacks the- 
oretical motivation. The pseudo-isothermal function 



p(r) = 



Po 



1 + {r/r c f 



(4) 



contains a central DM space density of po and scale length of 
r c . In the limit of r c significantly larger than the datapoints, a 
power-law function with a = mimics this function. How- 
ever, the added flexibility of the scale length term can dimin- 
ish the capability to discriminate between cores and cusps. 
The best-fit, five-parameter model with the SPS penalty has a 
minimum at T*„ R =1.37, i =66.4°, j3 z =0.390, p =0.198 
M Q pc~ 3 , and r c =1.0 kpc at xt=s8 =78.3. The steep, 
power-law model is still preferred, but the cored model is sta- 
tistically viable. Data at larger radii are necessary to better 
test for the presence of a large core. 

Next, we test a six parameter model with the NFW form. 
A minimization of the full NFW function requires a large 
scale radius, r s , for a quality fit. Only models with r s > 10 
kpc fit well where the power-law approximation becomes 
highly precise. This results in an enormous DM virial mass 
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(~ 2 x 10 n MQ), but data at larger radii are necessary to make 
a reliable estimate of the virial mass. 

5. MODELS FROM GAS KINEMATICS 

The [Oil] data are fit with a tilted ring (TR) dRogstad et all 
11974 119761) and harmonic decomposition (HD) algorithm to 
determine a rotation curve assuming an infinitely thin geom- 
et ry for the g a s. Th e code is the same as used and described 
in lFathi et al.l ([2005). The HD model and terms are shown in 
Figure [14] as is the TR model. Driven by the same complex 
structures and kinematic twists as discussed by SBLB03 for 
the Ha in this galaxy, we have make harmonic fits through 
the to = 3 terms. SBLB03 only present circular and radial 
terms, however. Our TR fit is allowed a position angle that 
varies with radius which can also explain the kinematic twist. 
The [Oil] rotation curves are given in Table [3] and shown in 
Figure Q3] The asymmetric drift correction for the ionized 
gas in NGC 2976 has been calculated in SBLB03, found to 
be small, and not used in their analysis because of substantial 
uncertainties in its exact value. Similarly, we do not apply 
an asymmetric drift correction to our gas rotation curve fits. 
There is consistency between the shapes of our [Oil] rotation 
curve and the Ha rotation curve of SBLB03. The irregular 
structure in the curve at r= 30" and r= 60" is found in both 
datasets, particularly in our tilted ring fit. 

We next fit a velocity power law of the form v c i rCj DM oc T ^ 
added in quadrature to the stellar and HI component circular 
velocities. The circular velocity of a power law density pro- 
file is a power law with a different index. The well-known 
relation for the density and circular velocity indices for power 
laws of a = 2 X (1 — j3) (e.g. SBLB03 Appendix B) is used. 
A variety of T* values are tested as detailed in Table|5J some 
fixed and some fit live. We also refit the SBLB03 data of their 
Table 3 in the same manner, taking their mass model rotation 
curves and trying their maximal disk value of M*/L#- = 0.19. 
The residuals from the best fit were used to estimate the un- 
certainty in the rotation curve; these systematic uncertainties 
are larger than the statistical errors and included in the error 
determination of a. 

The estimates of a from the gas kinematics are presented in 
Table|5]for a variety of datasets and assumptions on T*. The 
HD fits to the [Oil] data require a cored DM halo, regardless 
of the T* assumptions. The [Oil] TR fits require DM slopes 
that are steeper than the HD fits at 3a significance, but still 
deviating from NFW expectations at 1. 5a significance. The 
harmonic decomposition can fit for radial infall or outflow, but 
some other motions may be degenerat e with rotation such as 
the m otions due to a bar-like potential ( Snek kens & Sellwoodl 
20071). A bar is expected to show power in the third order 
terms, and our m = 3 sine term does show some power at 
large radius. Finally, we fit the harmonic decomposition data 
of SBLB03 in the same manner. From their fits with a range 
in M/L chosen to represent maximal and submaximal disks, 
SBLS03 reach ranges of 0.01 < a <0.17. We find agreement 
with their determinations. We find a larger, but still deviat- 
ing from a cusp by 5a significance, DM halo slope when we 
constrain the disk to have no mass. We conclude that a mass 
model based on the gas data with a single position angle and 
harmonic terms favors a cored DM halo, but that this result 
is not robust against an equally viable model that fits the data 
with a twist in the position angle. The mass profile conclu- 
sions drawn from the gas kinematics in this object are subject 
to severe dependencies in the modelling choices. 

A comparison between the gas velocity field in Figure [14] 



and the stellar line-of-sight velocity in Figure [6] shows that 
the zero velocity contours in both are twisted in the same way 
and with similar magnitude along the minor axis. This may 
be an important clue to the cause of the non-regular motions 
as collisional processes are expected to not twist the stellar 
zero velocity contour. One explanation may be that both the 
stars and gas have their angular momentum vectors perturbed 
at small radius, perhaps by bending in a weak bar potential. 
In this case, the stellar models are more immune to a warp 
as the Jeans models compare to V rms instead of simply Vj os , 
and we measure a stellar velocity dispersion that is larger than 
line-of-sight rotation out to R maj ~ 60". With the current 
level of stellar, observational errors, the axisymmetric models 
we present are statistically sufficient to describe the galaxy. 
Better data, taken in the future, may merit analysis with non- 
axisymmetric orbit-based models. 

6. CONCLUSIONS 

We present two-dimensional maps of stellar and gaseous 
kinematics in the late-type dwarf galaxy NGC 2976. The- 
oretical models of collisionless, cold dark matter predict a 
cuspy dark matter halo to exist in low mass halos when gravity 
alone is simulated. Baryonic feedback processes have been 
proposed as a mechanism to create a cored halo, which ob- 
servations such as ours may constrain. The leverage of the 
stellar kinematics as a collisionless tracer is a major advan- 
tage to our work. We fit the stellar kinematics with an ax- 
isymmetric, semi-isotropic Jeans model to measure the DM 
profile and constrain T*. The Jeans model permits both a 
cored or cuspy halo with a mild preference for a cored, high 
We next fit a suite of stellar population histories to an 
optical spectrum with a broader bandpass. We find a limit of 
T* fl=l.l±0.8 driven primarily by uncertainty in the initial 
mass function. This limit in combination with the kinematic 
data provides a much tighter certainty on the DM profile. The 
combined fit suggests a DM cusp (a = 0.90 ± 0.15 at la) 
and excludes a DM core at 2a significance. NGC 2976 is 
dark matter dominated everywhere outside of the nuclear star 
cluster and requires no history of baryonic feedback or non- 
standard particle properties to explain the dark matter halo 
profile. The gaseous kinematics, in concordance with ear- 
lier work, imply a cored dark matter halo when modeled with 
a harmonic decomposition method and a single position an- 
gle. A tilted-ring fit with a variable position angle instead 
excludes the cored model and is compatible with the mod- 
els fit by the stellar kinematics. A larger sample with stellar 
kinematics is motivated to compare with the vast literature on 
gaseous kinematics and reassess the quantity and distribution 
of dark matter and possible correlations with central density 
slope in late-type dwarf galaxies. This work comes with two 
primary caveats that can be improved upon. First, our results 
are strongest for a DM density power-law approximation, and 
a DM density model with a large core under a different param- 
eterization ca nnot be strictly excluded with the current data as 
discussed in i )4.5l Data at larger radii, such as >2 kpc, will 
close this uncertainty but require a deeper flux limit. Second, 
a crucial barrier to achieving lower statistical errors in this 
work has been the large instrumental resolution relative to the 
intrinsic dispersion. In future works, we will use the newly 
available Visible Integral field Replicable Unit Spectrograph 
Wendelstein (VIRUS-W) (iFabricius et al J 120081) at a resolu- 
tion of R = 6800 to circumvent this limitation. 
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Figure 2. The simultaneous wavelength and instrumental resolution fit for a particular VIRUS-P fiber. The solutions are parameterized by a fourth order 
polynomial in wavelength and a third order polynomial in spectral resolution. The fits are made by matching a higher resolution solar spectrum to twilight flat 
frames. 
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Figure 3. The instrumental resolution in each kinematic extraction bin. 
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Figure 4. The extracted velocity dispersions in simulated data. The instrumental resolution was set to 50 km s _1 , and the color is coded by stellar velocity 
dispersion. Top: The simulated data came from a subset of the LOSVD template set. No systematics are seen down to very low S/N and dispersions. Bottom: 
The extracted velocity dispersions in simulated data. The simulated data come from stars in the ELODIE dataset which are not contained in our LOSVD template 
set. This tests for template mismatch. No systematics are seen down to very low S/N and dispersions. 
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Figure 5. Spectra and fits in several bins. The chosen bins (56, 38, 19, and 5), represent a range of S/N (57, 57, 53, and 26) starting from the top left going 
clockwise. The error normalized residuals are shown in the bottom panel. 
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Figure 6. The velocity maps for NGC 2976. Top left: The stellar line-of-sight velocity map. Bottom left: The line-of-sight velocity error map. Top right: 
The stellar velocity dispersion map. The data are consistent with a flat or only modestly sloped velocity dispersion profile. Bottom right: The stellar velocity 
dispersion error map. 
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Figure 7. The rms velocity as measured in various spectral windows. The values are color coded by their distance along the major axis. The first, second, 
and third labels refer to the 3900A< A <4100A, 4100A< A <434()A, and joint windows respectively. The estimates generally agree, especially for the high 
S/N datapoints which drive the modeling. We use the 4100A< A <434()A kinematics as the preferred value due to the lower errors and the more homogenous 
instrumental resolution compared to the joint range. 




Figure 8. Data and fits to synthetic stellar population models. The spectra have been corrected for their radial velocity, broadened appropriately, and stacked. 
The SINGS spectrum has poorer resolution than the SPS models, so the SPS models are broadened to match. The average S/N per pixel over the displayed data 
is 60. Regions with emission lines have been masked in the fit and shaded in grey hatches. The best two-component model is shown with confidence intervals. 
Our VIRUS-P data agree well in absolute and relative flux calibration with the SINGS spectrum. The SPS fit is made to the SINGS spectrum, preferred by its 
wider bandpass and better coverage of the 4000A break. The displayed best-fit model proscribes T ti n = 0.90 with luminosity weighted fractions of 0.91:0.09 
for a Z = 0.2Zq, E s (B-V)=0, rg fr = 5 Gyr, 12 Gyr age model and a Z = 0.2Zq, E s (B-V) = 0.60, instantaneous burst, 25 Myr age model respectively. 
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Figure 9. The stellar population fits to the optical and NIR d ata. The fits pr edict mass-to-light ratios that are consistent with, but more loo sely cons t rained , than 
the fit to the SINGS spectrum. The mass modelling of de Blok et al. (2008) uses the 3.6/im photometry and a color relation derived from Oh et al. (200^) for a 
value of T*. 3.6^771 = 0.66 and claims that NGC 2976 is baryon-dominated. Note that this relation was based on stellar populations aged to 12 Gyr and does not 
match the observed K-[3.6 (im] color of NGC 2976. Our models, specifically fit to the spectrophotometry of NGC 2976, disagree. We show the best-fit model 
and confidence interval with various two-component SPS models. The best-fit model shown proscribes T* r = 0.64 and T tj 3.6 Mm = 0.18 with luminosity 
weighted fractions of 0.67:0.33 for a Z = O.2Z0, E S (B-V)=0.40, instantaneous burst, 100 Myr age model and aZ = 0.2Zq, E s (B-V)=0, instantaneous burst, 
2.5 Gyr age model, respectively. 
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Figure 10. The sector photometry of NGC 2976, folded into one quadrant, with the Multi-Gaussian Expansion (MGE) fit also shown. The values in the upper 
right of the profile plot indicate the angle from the major axis, in degrees, for each displayed cut. Left: The HI distribution is strongly asymmetric and requires 
some negative mass components in the model's central regions for a quality fit. To maintain a physical density, we have instead made a poorer fit with only 
positive components. However, the total HI mass is significantly smaller than the stellar and DM mass, and its inclusion or exclusion does not strongly alter the 
DM halo fits. Right: The R-band data and fit. The nuclear stellar cluster and the break between inner and outer disk near r~70" on the major axis and r~25" on 
the minor axis are captured in the fit. The two star forming regions near r~70" on the major axis remain as residuals since they are asymmetric. 
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Figure 11. Left: V rms . Right: V rms error for the data and residual/error for the models. Row 1: Data. Row 2: The model at the least-squares minimum fit 
made to the kinematic data alone. This model disagrees with the SPS determination on T, h. The model is characterized by xt = s7 = 77.0. Row 3: The model 
at the least-squares minimum fit to the combined SPS and kinematic data. Including the SPS penalty, the fit is characterized by xt—ss = 77.1. Row 4: The model 
second-moment velocity map for a fit fixed toT t ^ = l.l,a = 0. 1, and the remaining parameters freely fit. This model represents the class of DM dominated 
and cored models that are excluded by our data. The penalty in \' 2 comes from model velocities that rise more slowly than the data along the major axis and 
excess model velocities at minor axis offsets. This fit is characterized by X^=89 = 93.9. 
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Figure 12. The enclosed mass profiles and confidence intervals. Top: Using the kinematic data alone. The Y*^ fit here exceeds that fit through stellar population 
synthesis analysis. With these poor constraints, NGC 2976 may be DM or baryon dominated to the largest radius measured. The median models have similar 
amounts of enclosed DM and stellar mass up to 1 kpc where the DM begins to dominate. Bottom: Using the combined SPS and kinematic data. The DM halo 
dominates at least to r > 200 pc and perhaps everywhere. 
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Table 1 

Kinematic template stars 



HYPERLEDA # 


HD# 


Type 


[Fe/H] 


Average weight 


Intrinsic a (km s x ) 


00672 


HD181214 


F8m 


-0.01 


0.086 


30.8 


01269 


HD058923 


Fom 


0.21 


0.185 





01286 


HD062509 


KOIIIb 


0.02 


0.006 


11.1 


01298 


HD068017 


G4V 


-0.41 


0.067 


10.6 


01322 


HD073667 


K1V 


-0.55 


0.032 


10.9 


01359 


HD088609 


G5IIIwe 


-2.67 


0.197 


12.3 


01366 


HD089744 


F7V 


0.11 


0.008 


15.3 


01733 


HD169985 


GOIII+ 


0.34 


0.138 


12.3 


01828 


HD193793 


WC+ 


0.77 


0.280 






Table 2 

Stellar kinematic measurements in NGC 2976* 





Bin #' 


Major axis 
distance 
(") 


Minor axis 
distance 
(") 


^los,* 
(kms- 1 ) 


Vlos,*, 1 O 

uncertainty 
(kms" 1 ) 


cr * 
(kms- 1 ) 


cr*, 1 a 
uncertainty 
(kms- 1 ) 


1 


-48.1 


-62.5 


-34.53 


26.29 


10.00 


49.63 


2 


-57.7 


-49.7 


-55.44 


15.52 


10.00 


34.29 


3 


-67.3 


-37.0 


-34.53 


18.81 


48.27 


24.73 


4 


-76.9 


-24.2 


-45.49 


3.90 


10.00 


18.03 


5 


-86.5 


-11.4 


-51.80 


5.01 


10.00 


19.68 


6 


-96.2 


1.3 


-62.00 


4.67 


10.00 


19.78 


7 


-105.8 


14.1 


-65.17 


7.47 


10.00 


23.18 


8 


-35.3 


-52.9 


27.51 


9.08 


10.00 


26.72 


9 


-44.9 


-40.1 


-1 1.36 


5.43 


24.54 


12.43 


10 


-54.5 


-27.3 


-14.32 


4.37 


30.58 


7.89 


11 


-64.2 


-14.6 


-32.05 


3.31 


29.10 


6.55 


12 


-73.8 


-1.8 


-43.70 


3.18 


46.88 


5.49 


13 


-83.4 


11.0 


-42.71 


4.22 


22.76 


10.50 


14 


-93.0 


23.8 


-45.30 


10.15 


56.69 


17.22 


1j 


o o 
-3.3 


-05.0 


OA 1 o 

-2o.l5 


Ij.JO 


10.00 


OO 1 A 

51. 14 


16 


-12.9 


-56.0 


-6.00 


9.61 


39.36 


1 1.88 


1 / 


oo c 
-22. J 


A Q O 

-43.2 


-1.24 


C 1 A 
0.14 


1 A oo 
lO. 55 


12.89 


i o 
15 


oo 1 
-32.1 


on c 
-30. J 


-9.43 


O TO 

2. 13 


OA n 1 
20.91 


7.21 


19 


-41.8 


-17.7 


-18.39 


1.71 


30.71 


5.34 


20 


-51.4 


-4.9 


-26.85 


1 /lO 

1.43 


31.07 


4.85 


21 


-61.0 


7.9 


-29.12 


1.50 


34.13 


4.89 


22 


-70.6 


20.6 


-25.46 


3.28 


39.24 


5.65 


23 


-80.3 


33.4 


-39.14 


9.02 


43.85 


13.96 


24 


9.5 


-59.1 


-1 1.57 


6.45 


15.23 


13.59 


25 


-0.1 


-46.4 


-12.83 


4.69 


31.38 


7.91 


26 


-9.7 


-33.6 


1.99 


2.61 


22.52 


7.57 


27 


-19.4 


-20.8 


-8.08 


1.64 


32.75 


4.94 


28 


-29.0 


-8.0 


-12.27 


1.40 


32.78 


4.62 


29 


-38.6 


4.7 


-14.64 


1.35 


34.33 


4.71 


30 


-48.3 


17.5 


-12.42 


1.73 


32.77 


5.03 


31 


-57.9 


30.3 


-9.70 


3.21 


45.44 


5.29 


32 


-67.5 


43.1 


-7.36 


8.60 


48.44 


14.59 


33 


31.9 


-62.3 


-11.61 


5.59 


35.13 


10.25 


34 


22.3 


-49.5 


-1.60 


3.29 


33.96 


7.72 


35 


12.7 


-36.7 


1.96 


2.54 


22.54 


7.65 


36 


3.0 


-24.0 


1.11 


1.85 


26.62 


5.67 


37 


-6.6 


-11.2 


-6.59 


1.73 


25.43 


5.62 


38 


-16.2 


1.6 


-7.58 


1.54 


27.82 


5.04 


39 


-25.9 


14.4 


-4.36 


1.32 


32.51 


4.59 


40 


-35.5 


27.2 


-7.92 


1.84 


32.87 


4.89 


41 


-45.1 


39.9 


-10.69 


3.90 


25.57 


8.49 


42 


-54.7 


52.7 


-26.01 


8.98 


38.92 


16.60 


43 


44.7 


-52.7 


9.80 


5.35 


41.22 


9.30 


44 


35.1 


-39.9 


7.40 


3.25 


24.71 


7.73 


45 


25.4 


-27.1 


3.85 


1.69 


23.17 


6.58 


46 


15.8 


-14.3 


5.60 


1.86 


31.52 


4.77 


47 


6.2 


-1.5 


6.45 


1.43 


30.24 


4.65 


48 


-3.5 


11.2 


5.29 


1.50 


28.17 


4.99 


49 


-13.1 


24.0 


7.00 


1.74 


28.55 


5.06 


50 


-22.7 


36.8 


5.31 


1.93 


29.17 


6.11 


51 


-32.4 


49.6 


-2.15 


5.02 


18.75 


12.36 


52 


-42.0 


62.4 


-7.00 


19.10 


45.28 


28.29 


53 


57.5 


-43.0 


11.07 


3.43 


38.57 


6.14 


54 


47.8 


-30.2 


15.94 


2.74 


31.43 


6.30 


55 


38.2 


-17.5 


13.44 


1.52 


25.29 


5.70 


56 


28.6 


-4.7 


16.54 


1.54 


34.99 


4.28 
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Table 2 — Continued 



Bin #t 


Major axis 


Minor axis 


v LOS,* 


V; 1 rj 

V fOS,*' 1 u 




(7* , 1 (T 




Histanre 


Li.±L> I til 


^Km S ) 


linrertaintv 




nnrprtnintv 




I ) 




I Kill S ) 


^KIll S ) 


^KITl S ) 


57 


18.9 


8.1 


20.57 


1.43 


26.04 


5.32 


58 


9.3 


20.9 


14.82 


1.53 


25.41 


5.78 


59 


-0.3 


33.7 


7.83 


2.53 


29.96 


5.51 


60 


-10.0 


46.4 


-1.65 


4.03 


34.12 


6.17 


61 


-19.6 


59.2 


34.47 


7.82 


27.31 


17.13 


62 


-29.2 


72.0 


43.28 


66.67 


10.00 


44.09 


63 


70.2 


-33.4 


25.36 


3.12 


27.98 


6.72 


64 


60.6 


-20.6 


25.45 


2.55 


42.12 


4.86 


65 


51.0 


-7.8 


22.35 


1.68 


36.79 


4.00 


66 


41.3 


5.0 


19.14 


1.46 


36.58 


4.02 


67 


31.7 


17.7 


18.32 


1.47 


33.43 


4.37 


68 


22.1 


30.5 


9.92 


1.62 


34.08 


4.60 


69 


12.4 


43.3 


4.18 


3.05 


31.55 


5.96 


70 


2.8 


56.1 


-5.93 


6.36 


10.00 


18.74 


71 


83.0 


-23.8 


44.14 


3.44 


39.98 


6.14 


72 


73.4 


-11.0 


37.15 


2.95 


10.00 


14.62 


73 


63.7 


1.8 


32.91 


2.67 


36.58 


6.40 


74 


54.1 


14.6 


32.88 


1.79 


27.56 


5.30 


75 


44.5 


27.4 


20.85 


2.07 


32.52 


4.91 


76 


34.8 


40.2 


7.99 


2.92 


37.76 


5.12 


77 


25.2 


53.0 


5.11 


4.08 


10.00 


15.89 


78 


15.6 


65.7 


-7.83 


6.33 


23.54 


11.61 


79 


95.8 


-14.1 


44.77 


3.09 


24.18 


7.57 


80 


86.2 


-1.3 


44.70 


2.81 


26.93 


6.62 


81 


76.5 


11.4 


37.63 


3.08 


39.61 


4.68 


82 


66.9 


24.2 


31.54 


2.88 


45.35 


4.49 


83 


57.3 


37.0 


31.35 


4.23 


38.50 


6.47 


84 


47.6 


49.8 


19.65 


7.23 


25.00 


13.79 


85 


38.0 


62.6 


-5.33 


11.58 


37.91 


17.95 


86 


108.6 


-4.5 


38.98 


7.02 


27.43 


11.96 


87 


98.9 


8.3 


43.17 


7.08 


44.62 


9.69 


88 


89.3 


21.1 


22.96 


5.53 


40.10 


8.91 


89 


79.7 


33.9 


40.93 


12.32 


48.21 


16.14 


90 


70.0 


46.7 


2.96 


8.51 


39.16 


17.83 


91 


60.4 


59.4 


-16.52 


8.01 


35.10 


19.13 


92 


50.7 


72.2 


20.40 


40.00 


25.82 


24.16 



* Column 4 contains the systemic velocity which we measure as a weighted average to be 4.60 km s~ 1 . The quantity fit in the Jeans models is the velocity second-moment (V rms ), 
or the sum in quadrature of columns 4 and 6 with the systematic velocity removed. The V rTTls error is formed by the propagation of the listed terms as well as an estimated 0.5 km 
s - 1 uncertainty on the average 50 km s — 1 instrumental resolution. 

^ The central position of each bin is given. The bins are 16" x 16"in size and rotated by -53° east of north. 



Table 3 

Rotation curve data 



Radius 


v c,TR,lOII]^ 


"c.bdjo;;]' 


v rad,[OII] 


Vc,** 


Vc,HI 


(arcsec) 


(kms' 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms- 1 ) 


(kms" 1 ) 


14.1 


15.6 ± 0.8 


10.6 ± 0.8 


5.2 ± 0.4 


9.5 


1.5 


19.1 


32.4 ± 0.5 


24.6 ± 1.0 


9.0 ± 0.6 


11.2 


2.0 


24.1 


39.4 ± 0.9 


31.4 ±2.5 


11.0 ±0.7 


13.0 


2.5 


29.1 


45.7 ± 0.8 


32.1 ± 1.6 


6.4 ± 1.5 


15.0 


3.0 


34.1 


40.0 ± 0.4 


32.4 ±5.9 


6.7 ± 1.6 


16.9 


3.5 


39.1 


43.5 ± 1.1 


39.7 ± 2.9 


7.4 ± 1.8 


18.8 


3.9 


44.1 


51.6 ±0.6 


44.7 ± 3.3 


2.0 ± 1.8 


20.5 


4.4 


49.1 


55.9 ± 0.7 


47.7 ± 3.3 


3.8 ± 1.3 


22.1 


4.9 


54.1 


59.7 ± 1.0 


51.7 ±0.8 


3.3 ± 1.0 


23.5 


5.3 


59.1 


56.2 ± 0.5 


51.0 ±2.8 


2.5 ± 5.5 


24.8 


5.7 


64.1 


57.3 ±0.2 


53.7 ±0.4 


-0.3 ± 0.9 


25.9 


6.1 


69.1 


63.4 ± 0.3 


61.8 ±0.8 


-2.2 ± 2.1 


26.8 


6.5 


74.1 


66.7 ± 0.4 


64.4 ± 0.8 


0.1 ± 1.1 


27.6 


6.9 


79.1 


69.9 ± 0.2 


72.4 ±4.8 


-4.0 ± 2.3 


28.2 


7.2 


84.1 


75.6 ± 0.5 


73.8 ± 0.9 


-3.5 ± 2.7 


28.7 


7.6 


89.1 


73.7 ±0.5 


75.5 ± 3.2 


-6.7 ± 5.3 


29.0 


7.9 


94.1 


71.3 ±0.6 


68.9 ± 1.6 


-4.2 ± 2.1 


29.2 


8.2 


99.1 


79.9 ± 0.5 


78.1 ± 13.8 


-17.2 ± 10.0 


29.4 


8.5 


104.1 


80.5 ± 0.9 


82.5 ± 8.4 


-14.8 ±7.7 


29.4 


8.7 


109.1 


74.5 ± 2.1 


69.0 ± 10.9 


-10.8 ± 3.8 


29.3 


9.0 



' The listed errors are statistical. The rms velocities from the best-fit TR and HD models are 4.5 km s _1 and 6.1 km s _1 for Y*,i? — 0, 4.4 km s _1 and 6.4 km s _1 for 
T„ H = 1.1, and 4.5 km s -1 and 6.0 kms -1 for a freely fit T„ iR , respectively. 
* Assuming T*,j? = 1. This column scales as oc (T.j) 1 ' 2 
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Figure 13. The A\ 2 distributions for combinations of the five parameters being minimized in the Jeans modelling. The inner three contours correspond to 
Ax 2 = 2.3, 6.2, and 12.9 or 1, 2, and 3 a significance. The additional contours increase each by a factor of two. Cross symbols mark the formal minima. Left: 
The confidence intervals from the kinematic data alone. Right: The confidence intervals when the kinematic and SPS data are combined. 




Figure 14. The data and fits to the [Oil] velocity field representing the two ways that the kinematic twist can be recovered. Left: The fit with a tilted ring (TR) 
model and a variable position angle. Right: The fit with a harmonic decomposition (HD) model through three orders. Signficant radial velocity is observed, 
particularly in the galaxy's center. The sign of the radial velocity is uncertainty as we cannot identify the galaxy's near side. If we speculate that the galaxy's near 
side is in the SW, the gas near the center is outflowing. The higher-order terms are noisy, but similar rotational and radial velocity fits are made when the higher 
order terms are not fit. 



Table 4 

Multi-Gaussian Expansion terms for NGC 2976 



Component 


Index 
(k) 


(M Q pc" 2 ) 


(arcsec) 


* 


Stellar* 


1 


2298.3 


0.28 


1.00 




2 


68.9 


0.99 


1.00 




3 


27.8 


6.36 


1.00 




4 


117.8 


49.89 


0.48 




5 


19.6 


108.26 


0.54 


HI 


1 

2 


8.0 
1.3 


84.80 
291.13 


0.70 
0.75 



R-band and assuming Y* = l for this listing. 
Observed axial ratio. 
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Figure 15. Comparison between the circular velocity determinations from stellar and gas kine mat ics. The solid line in each panel shows the circular velocity 
for the mass model favored by the stellar kinematics with a SPS penalty (right panel of Figure [T3V The circular velocity as determined by the gas kinematics 
is shown with datapoints and the dotted line as a power-law fit. The 68% confidence interval is indicated by an upper and lower dotted line. The fits to the gas 
data use fixed Y„ values per Table[5] Left: The [Oil] data and fit for the tilted -ring model with a varying position angle. Middle: The [Oil] data and fit for the 
harmonic decomposition model through third order. The fit is steeper than the stellar based model and consistent with a cored DM halo. Right: The Ha data of 
SBLB03. These data are fit by removing a radial velocity component at small radius and similar to our [Oil] harmonic decomposition model. 

Table 5 

Dark matter density index constraints from gaseous kinematics t 



Assumption 


This dataset's 


This dataset's 


SBLB03 dataset 




tilted ring fit 


harmonic decomposition 




r» =o 


0.80 ± 0.08 


0.31 ± 0.13 


0.43 ± 0.06 


T, fixedtt 


0.84 ±0.10 


0.14 ±0.17 


0.16 ± 0.08 


freely fit T, 


0.81 ± 2.20 


0.31 ±0.98 


0.12 ±0.44 



The confidence ranges are based on rms circular velocity uncertainties in order to capture both statistical and systematic errors. 
' Assuming T * i ^ — 1.1 for our data and T* k =0.19 for the SBLB03 data. 



